Bayesian analysis of inflation: Parameter estimation for single field models 
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Future astro-physical datasets promise to strengthen constraints on models of inflation, and ex- 
tracting these constraints requires methods and tools commensurate with the quality of the data. 
In this paper we describe ModeCode, a new, publicly available code that computes the primordial 
scalar and tensor power spectra for single field inflationary models. ModeCode solves the inflation- 
ary mode equations numerically, avoiding the slow roll approximation. It is interfaced with CAMB 
and CosmoMC to compute cosmic microwave background angular power spectra and perform like- 
lihood analysis and parameter estimation. ModeCode is easily extendable to additional models 
of inflation, and future updates will include Bayesian model comparison. Errors from ModeCode 
contribute negligibly to the error budget for analyses of data from Planck or other next generation 
experiments. We constrain representative single field models ((p™ with n — 2/3, 1, 2, and 4, natural 
inflation, and "hilltop" inflation) using current data, and provide forecasts for Planck. From current 
data, we obtain weak but nontrivial limits on the post-inflationary physics, which is a significant 
source of uncertainty in the predictions of inflationary models, while we find that Planck will dra- 
matically improve these constraints. In particular, Planck will link the inflationary dynamics with 
the post-inflationary growth of the horizon, and thus begin to probe the "primordial dark ages" 
between TeV and GUT scale energies. 



I. INTRODUCTION 



The last two decades have seen a sequence of break- 
throughs in the understanding of the physical universe. 
The detection of cosmic microwave background (CMB) 
anisotropies by COBE in 1992 [l[, evidence for dark en- 
ergy in the distance-luminosity relationship for Type la 
supernovae in 1998 @, Q, and the sequence of WMAP 
data releases [4H7[, among others, mark turning points 
in our ability to constrain — and falsify — specific cos- 
mological models. These advances begin to fulfill the 
long-standing promise that astrophysical data can di- 
rectly probe the first moments after the Big Bang, while 
simultaneously constraining models of ultra-high energy 
physics. Experiments underway and now being devel- 
oped guarantee that this progress will continue well into 
the future. In particular, the Planck satellite [8| has com- 
pleted a full survey of the sky, and this data should dra- 
matically improve the constraints on the free parameters 
in the so-called concordance cosmology. 

Given the quality of present-day data, the primordial 
perturbations are fully described by two numbers (e.g. 
Ref. Q): the amplitude A B and tilt n s of the power 
spectrum of density (scalar) perturbations. These quan- 
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tities form the basis of an inflationary sector in the con- 
cordance model, if we postulate that their values can 
be traced back to a phase of primordial inflation. Con- 
versely, a complete inflationary model must account for 
these numbers. As the data improve, this set of param- 
eters can easily expand to include the properties of any 
primordial gravitational wave background or departures 
from Gaussianity. Given the tight agreement between the 
concordance model and current data, the impact of any 
new parameters is necessarily subleading. Consequently, 
additional parameters needed to describe the primordial 
perturbations can be regarded as fingerprints of specific 
inflationary models, in that most of these quantities will 
be vanishingly small in most inflationary models (see e.g. 
Ref. [3). 

Physically, inflation is characterized by a period of ac- 
celerated expansion in the early universe, and an infla- 
tionary model is defined by the mechanism that drives 
this accelerated expansion. Simple models of inflation 
can usually be described by the kinetic term and poten- 
tial of a single scalar degree of freedom (the inflaton), 
along with this field's coupling to gravity. In this paper 
we focus on models where the field is minimally coupled 
to gravity and has a canonical kinetic term, but will relax 
these restrictions in future work. 

The simplest approach to constraining inflation is to 
specify the primordial perturbations in terms of the em- 
pirical parameters 
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where a s is the running of the scalar spectral index, r the 
tensor-to-scalar ratio, n t the tensor spectral index, and 
/ n i parametrizes non-Gaussianity. Further, 6 cmp may be 
extended indefinitely to include higher-order terms in the 
expansion of the scalar and tensor power spectra and 
various deviations from Gaussianity. These additional 
parameters will generally be more difficult to measure in 
any given dataset than the basic empirical parameters A s 
and n s . In all cases, inflationary model predictions for 
the values of the empirical parameters are to be compared 
with the measured values of these parameters. 

A second approach treats the determination of the in- 
flationary mechanism as an inverse problem, and thus 
reconstructs the inflaton potential from the data [11H15] . 
For instance, slow roll reconstruction [l6Hl9( uses a sys- 
tematic expansion based on the slow roll hierarchy [20|— 
|25| together with consistency conditions on the duration 
of inflation, providing a minimally parametric approach 
to the inverse problem. Variations to this scheme have 
been widely discussed in the literature J26l433l |. The in- 
verse problem does not have a unique solution, and typ- 
ically encodes basic assumptions regarding the general 
class of inflationary models under consideration 34 1 . 

Here we adopt a third approach, in which we assume 
that inflation was driven by a specific model (i.e. form of 
the potential) with one or more free parameters, 



v = {Vi,V 2 ,...,V n ;d Ka }, 



(2) 



where the Vi parametrize the potential, while 0rh 
parametrizes the post-inflationary reheating phase, as 
we discuss later. These parameters replace the empiri- 
cal inflationary sector, omp . We estimate their values 
alongside other cosmological parameters, typically using 
Markov Chain Monte Carlo (MCMC) analysis jaJJEJI. 
Following parameter estimation for several different in- 
flationary models, Bayesian model selection p3 - |47j tech- 
niques will allow us to compare the fidelity with which 
these different models account for observations of the sky. 

Working directly with inflationary models, we can 
compute the power spectrum numerically without re- 
course to the slow roll approximation. This removes a 
lingering uncertainty from the analysis of the inflationary 
parameter space, and allows us to constrain inflationary 
models with complicated spectra not well described by 
9 cmp . While the exact computation of the power spec- 
trum is numerically straightforward and has been used 
extensively for slow roll-violating potentials previously 
[3ll . |32| . |48to9| , such precision is only now becoming nec- 
essary with the arrival of higher quality data from Planck. 

Even when the slow roll approximation is accurate and 
the spectra simple, there are advantages to working with 
By, rather than 9 cmp . Empirical quantities such as A s , 
n s , r, etc. are computable functions of the free param- 
eters of any specific inflationary model, and thus con- 
straints on C mp can, in principle, be mapped into con- 
straints on Oy . However, we will see that this approach 
is of limited use in practice, while the opposite approach 
of inferring constraints on cmp from constraints on Oy 



(given the underlying inflationary prior) is generally re- 
liable. 

Finally, By explicitly accounts for the theoretical un- 
certainty in inflationary predictions induced by the un- 
known thermalization history of the post-inflationary 
universe [17], UM, Io0l - l63l ] , in combination with the scale- 
dependence of the spectral index. This uncertainty is 
significant even when the running a s = dn s /dlnk is not 
itself observable. This apparent paradox arises because 
when the running is included, cmp requires at least four 
free parameters to fix the primordial spectra. Conversely, 
9y may have only a single free parameter for the infla- 
tionary potential and another for the reheating physics, 
and these numbers can thus be determined with more 
precision than those in e mp [64[ ■ 

This paper has three primary objectives. The first is to 
introduce ModeCode, a plug-in for CAMB and Cos- 
MOMC [H, Hlj]. ModeCode provides an efficient and 
robust numerical evaluation of the inflationary pertur- 
bation spectrum, and allows the free parameters in the 
potential to be estimated within an MCMC computation. 
Secondly, we use this code to generate constraints on rep- 
resentative single field inflationary models using current 
data, and give forecasts for the constraints that can be 
expected from Planck. Finally, this analysis underlines 
the importance of assumptions regarding reheating and 
thermalization in studies of inflation using data from the 
next generation of astrophysical datasets. We will see 
that current data put weak but nontrivial constraints 
on the reheating history given an explicit inflationary 
model. Further, our forecasts suggest that Planck will 
link the post-inflationary history with the inflationary 
epoch, with significant implications for theories of parti- 
cle physics between Grand Unified Theory (GUT) scales 
(~ 10 15 GcV) and TcV scales. 

ModeCode computes both the scalar and tensor per- 
turbation spectra, via the algorithm described in Ref. 
[50(1 . A number of common inflationary potentials are 
already included, and new models are straightforward to 
incorporate in ModeCode by supplying the potential 
and its derivatives. Since the CMB angular power spectra 
and likelihood are already expressed numerically, noth- 
ing is lost (other than a relatively small computational 
overhead) by solving directly for the mode amplitudes, 
rather than using an analytic approximation. Attention 
has been given to ensuring that the initial conditions 
for the background arc self-consistent and that the code 
"fails gracefully" for unphysical parameter combinations, 
so that such points are excluded from MCMC analyses. 
In addition, we accurately compute the endpoint of in- 
flation and the evolution of the comoving horizon size 
during inflation, so as to precisely match scales in the 
inflationary era with scales in the present-day universe. 

In a follow-up paper, we will interface ModeCode 
with MultiNest [431, allowing us to compute the 
Bayesian evidence for the models we analyze. In addi- 
tion, it is straightforward to extend ModeCode to sce- 
narios with nonminimal kinetic terms 15a. 15711 , multi- field 



models, or even non-Gaussianity |6£ 



II. METHODS 



background equations in terms of In a. Denoting d/dhia 
with a prime (') and recalling that H = dhia/dt by def- 
inition, we write the Einstein equation for H' and the 
Klein-Gordon equation for <j> as follows: 



A. Numerical solution 

In many circumstances, the slow roll approximation 
provides a sufficiently accurate description of the infla- 
tionary power spectra and has the advantage of express- 
ing them as functions of the potential and its derivatives. 
However, we want to avoid the slow roll approximation to 
maintain full generality and accuracy, and instead com- 
pute the initial curvature and tensor power spectra nu- 
merically given a specific inflaton potential V(<f>) |48| - |50| . 
This approach yields exact numerical results for arbitrary 
inflationary potentials, up to the intrinsic accuracy of 
first order gravitational perturbation theory. 1 

We begin by reviewing the formalism used for the nu- 
merical solution of the mode equations. We describe the 
scalar perturbations with the gauge invariant Mukhanov 
potential u j67|, |68[, which is related to the curvature 
perturbation 1Z: 



u = —zTZ . 



(3) 



where z = 4>/H, H is the Hubble parameter, and dots 
denote derivatives with respect to conformal time. The 
Fourier components Uk obey [69h71| 
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0. 



(4) 



where k is the modulus of the wavevector k. The power 
spectrum is defined in terms of the two point correlation 
function 



(ft k ft£,) 



k 3 



and is related to Uk and z via 



A 2 



*<*> - 2^ 



In terms of the empirical parameters A s , n s , and a s 



A%(k) = A s 



^pivot 



-l+^a s ln(fe/fcpi vot ) + - 



(6) 



(7) 



where fc p i vo t denotes the pivot scale at which the power 
spectrum is normalized. 

Equation (jj) depends on the background dynamics 
through z and its derivatives. Since the logarithm of the 
scale factor is a natural time coordinate for numerical so- 
lutions of the inflationary mode equations, we express the 



1 Second and higher order mode-mode couplings lead to both non- 
Gaussianity and loop corrections to the two-point functions. For 
the models discussed here, these corrections are very small. 



H' = - 



M 2 
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H + 



b') 2 H, 



1 dV 



= 0, 



(8) 
(9) 



where Mpi is the reduced Planck mass. Conveniently, 
our choice of independent variable gives </)' = z, and with 
the help of these background equations, the mode equa- 
tion ((3]) can be written 
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-H +1 






a 2 H 2 
H' 1 d 2 V 



2-4 — — 
H 4>' 



H H 2 d(j> 2 



u k = 0, (10) 



where the term in square brackets is z / (za 2 H 2 ) . 

We begin the integration of the background equations 
when the mode of interest is still deep inside the horizon 
(i.e. k 3> WOaH). We set the initial field velocity to its 
slow roll value, solving only the background equations, 
ensuring that the (small) initial transient in the velocity 
is damped away. When the mode is roughly l/100th of 
the horizon size we start to evolve the two orthogonal so- 
lutions that contribute to Uk , and read off the asymptotic 
value of \uk/z\ when the mode is far outside the horizon 
and frozen, as explained in Ref. J50(. 

The usual mode equation for tensor perturbations, 



t'k 



A 2 z (k)(27r) 3 S^(k~k') 1 (5) becomes 



§ + ' 



v k = , 



2 + ^ 
a 2 H 2 { H 



(11) 



v k = (12) 



after transforming the independent variable. The pri- 
mordial tensor power spectrum is 



A?(A) 



4 fc 3 



(13) 



and the asymptotic value of |ufe/a| is again taken from 
the numerical solutions. 



B. Matching and Theoretical Uncertainties 

Ideally, a complete model of the early universe would 
be predicted by a candidate theory of fundamental 
physics. In that case, we would know the mecha- 
nism by which energy is drained from the inflaton to 
(rc)thcrmalize the universe, as well as the equation of 



state and expansion rate of the primordial universe. Un- 
fortunately, inflationary model building is not mature 
enough for this to be done on a routine basis and, as a 
consequence, the unknown expansion history of the post- 
inflationary universe introduces a theoretical uncertainty 
into the predictions of inflationary models [64| . We must 
therefore introduce at least one phenomcnological param- 
eter #r,h to account for our ignorance of post-inflationary 
physics. For reasons explained below, we work with N, 
the number of e-folds of inflation between moment at 
which a specified mode leaves the horizon (k = aH) 
and the end of inflation, defined by the instant at which 
d 2 a(t)/dt 2 = 0. Further, given a specified pivot scale 
(fcpivot) we can then specify a corresponding iVpivot- For 
a given inflationary model. JVpivot is computed from the 
matching equation [6fJ, |72| . |73| . 

The matching depends on the growth of the horizon 
scale following inflation. This is a function of the detailed 
composition of the primordial universe, which is not well 
known, and the matching is thus intrinsically ambigu- 
ous. The two things we know with certainty are that 
the universe is not thermalizcd at the end of inflation, 
and that it must be thermalized by MeV scales, when 
primordial nucleosynthesis occurs. 2 Given that inflation 



can be a GUT scale phenomenon, the energy (r 



4/4 



may change by a factor of 10 18 between inflation and nu- 
cleosynthesis — a much larger factor than that between 
nucleosynthesis and the present day. 

A common assumption is that the universe is effec- 
tively matter-dominated (a(t) ex t 2 / 3 ) as the inflaton 
oscillates about the minimum of the potential, ther- 
malizes at some temperature Trh (with prh °c Tg H ), 
and remains radiation-dominated (a(t) ex t 1 ' 2 ) until the 
matter-radiation transition at z k 3f50 [7j. Even for 
this assumed history, the uncertainty in Trh allows a 
wide range in N. For example, a universe which is effec- 
tively matter-dominated between energies of 10 15 GcV 
to 10 3 GeV needs 9 e-folds less inflation after a given 
scale crosses the horizon than one which thermalizes at 
10 GeV. In many simple inflationary models, the run- 
ning a s is a few times 10 , so the resulting uncertainty 
in the scalar spectral index n s is ANa s ~ 5x 10~ 3 , 
of the same order as the statistical error expected from 
Planck p, |64(. However, far more extreme possibilities 
for the post-inflationary expansion rate exist, including 
kination (a(£) oc t 1 / 6 ) [74J], frustrated cosmic string net- 
works (a(t) oc t) [751 ] , or even a short burst of thermal in- 
flation [76| , and taking these scenarios into consideration 
greatly magnifies the uncertainty in inflationary predic- 
tions for a given model. 

The truly fundamental variable which specifics the por- 
tion of the inflaton potential being traversed as the pivot 



mode leaves the horizon is simply 



Npivot 



the value of 



2 Recent evidence points to the existence of a cosmological neu- 
trino background Q, which freezes out at temperatures slightly 
higher than those that apply during nucleosynthesis. 



the inflaton as the fc p i vo t leaves the horizon. However, 
</> increases in some potentials and decreases in others, 
while the mapping <f> — > 4> + <f>o produces a new poten- 
tial with identical inflationary dynamics, so the numerical 
value of </> is not informative on its own. Conversely, N{(j>) 
has the same interpretation in all inflationary models and 
is a monotonic (and usually simple) function of <f> fo r a 
given potential. For these reasons we take iV p i vot rather 
than </> p ivot as a free parameter. The coupling between 
-Npivot and the post-inflationary universe simply reflects 
the physical reality that the observed inflationary pertur- 
bation spectrum is a function of the post-inflationary ex- 
pansion history, but our choice ensures that the primor- 
dial perturbation spectrum is calculated solely in terms 
of parameters that describe the inflationary epoch which 
generated it. 3 

As we pointed out above, iVpivot depends on the de- 
tailed expansion history of the post-inflationary universe. 
However, there are a huge number of possible combi- 
nations of phases in the early universe, and these can 
have strongly degenerate predictions for iV p i vot . 4 One can 
imagine that the unknown expansion history is replicated 
by an effective barotropic fluid with equation of state 
w [63l . |64| , which is superficially equivalent to regarding 
IVpivot as a free parameter. Unfortunately, as explained 



in Ref. 6J], w is an ambiguous parameter. Given an ex- 



plicit inflationary potential, we can always determine the 
moment at which inflation ends, but we cannot compute 
w without specifying an energy scale at which the uni- 
verse has definitely thermalized: the numerical value of 
w is a function of this choice. Admittedly, A p i vot has an 
analogous dependence on the choice of fc p i vo t- However, 
assuming slow roll, iV p i VOt ~ log (/c p i VOt )/(l — e) so this 
dependence is usually transparent. Further, for a given 
combination of datasets, it is possible to determine an 
optimal choice of fc p i VOt (l7| . which is typically close to 
the geometric mean of the range of scales contributing to 
the dataset(s). 

We note that there is a correlation between the spectral 
amplitude at fc p i V ot (^4 S in the usual ACDM parameter 
set), and iV p i vot . As an example, consider m 2 (j) 2 infla- 
tion: to first order, <^ e nd and N((j>) do not depend on m 2 , 
while lowering m 2 lowers A s for fixed <j>. Consequently, 
V cn d will also decrease, ensuring that slightly less growth 
occurs between the end of inflation and some fixed ref- 
erence point, such as nucleosynthesis or recombination. 
Given the precision with which the spectral amplitude 
is now measured, this effect is small. Moreover, a simi- 



3 This approach works for scalar perturbations, since modes at 
astrophysical scales today remain outside the horizon until af- 
ter nucleosynthesis, at which point the thermal history of the 
universe is well constrained. However, primordial gravitational 
waves seen by direct detection experiments can enter the horizon 
during epochs for which the expansion history is not tightly con- 
strained, and for these we need the full transfer function l77tl78H . 

4 For example, a long matter-dominated phase leads to the same 
prediction for N p - lvot as a suitable combination of early radiation- 
domination and a short secondary period of inflation. 



lar degeneracy arises with w, since this parameter is also 
sensitive to changes in the energy density at the end of 
inflation, if all other parameters are held fixed. 

Theoretical considerations put very broad constraints 
on iVpivot- Firstly, in order to ensure that modes do ac- 
tually reenter the horizon, we need w > —1/3, so that 
p + 3p > and a < 0. This is not incompatible with a 
secondary period of inflation, but does require that the 
average expansion is not inflationary. Secondly, for a 
barotropic fluid, p > p is required in order to avoid a 
superluminal sound-speed, so w < 1. We are free to es- 
timate inflationary parameters for a narrower range of w 
or iVpivot, but it is important to recognize that doing so 
amounts to imposing a theoretical prior on the properties 
of the post-inflationary universe. 

In what follows, we constrain the free parameters of 
inflationary potentials for two different reheating sce- 
narios, general reheating (GRH) and instant reheating 
(IRH). In the latter case, we assume that the universe 
thermalizes instantaneously as inflation ends, and re- 
mains thcrmalizcd until matter-radiation equality In 
contrast, GRH assumes only that the universe is ther- 
malized by nucleosynthesis scales and that the aver- 
age expansion is no slower than radiation-dominated, or 
N < /V (IRH) 

iVnivot S JV pivot , 



v pivot 



the value computed assuming instanta- 
neous thermalization. This prior is analogous to stipu- 
lating that —1/3 < w < 1/3, which implicitly rules out a 
long kination-like phase. Lower values of w correspond to 
smaller N p - lvo t, which for the models considered here leads 
to an increasingly red-tilted spectrum. Current data are 
thus more sensitive to lower values of w for the models 
considered here. 



C. MCMC methodology 

Our Markov Chain Monte Carlo methodology is 
based on a modified version of CosmoMC [39j . using 
Metropolis-Hastings sampling for basic parameter esti- 
mation, and nested sampling via a MultiNest J47( plug- 
in for the calculation of Baycsian Evidence (to be pre- 
sented in a forthcoming publication). The free parame- 
ters in the potential V((f>) (plus the one reheating param- 
eter in the GRH case) are varied in the MCMC chains, 
along with the other cosmological parameters of the con- 
cordance model, and any nuisance parameters associated 
with the datasets. 

For all MCMC analyses of current CMB data, we run 
> 6 chains per model/data combination, requiring the 
Gelman- Rubin [79| criterion on the eigenvalues of the co- 
variance matrix to be R — 1 < 0.01 for convergence. For 
Planck simulation runs we use 4 chains per model which 
satisfy R — 1 ~ 0.1. While this convergence criterion is 
not as rigorous as that used in our main analysis, we ex- 



pect the resulting uncertainty in the estimated variances 
to be exceeded by the foreground-removal uncertainties 
[80J in the large-angle i?-mode constraint, which are not 
taken into account in these forecasts. 



D. Initial conditions and reheating 

We will now describe the implementation of the infla- 
tionary initial conditions and the reheating scenarios in 
ModeCode. 

Initial conditions: Inflationary potentials differ in their 
sensitivity to initial conditions (see Ref. [8l| and refer- 
ences within). Thus, automatically setting self-consistent 
initial conditions is a nontrivial issue. The initial value of 



</>' is set according to the slow roll equations (i.e. the in- 
flaton is assumed to be initially on the slow roll attractor 
solution). In addition, for a given set of potential param- 
eters, the algorithm must produce a starting field value 
(/>i n it which corresponds to a time well before the modes 
of interest leave the horizon. For the particular models 
we consider here, the code produces a first guess for <^>i n it 
from the field value needed to achieve iV((/) init ) = 70 in 
the slow roll approximation: 



N(<kmt) 



1 



M& 



V 



(14) 



where </> on d is the field value at the end of inflation and 
V<f, = dV/d(j). The code then iterates on this initial guess 
until a self-consistent value of </>i mt is found. For some 
combinations of potential parameters, it is possible that 
no set of self-consistent initial conditions exists. In such 
cases we reject the parameter combination in the MCMC 
analysis by assigning it a very small likelihood. 

Reheating: ModeCode evolves the inflationary back- 
ground solution through to the end of inflation, defined 
by d 2 a(t)/dt 2 = which corresponds to 



(if 1 



(15) 



where en is the first Hubble slow roll parameter. 6 This 
calculation yields the number of e-folds between the ini- 
tial conditions and the end of inflation. The matching 
equation then gives the scale factor at the end of infla- 
tion, a cn d- We connect a physical "pivot" wavenumber, 
fcpivot, to a particular epoch during inflation using 



-N„ 



l, . n . II . n o — iV pivot TT 

"pivot — "pivot-*-* pivot — a cnd c LL 



pivot J 



(16) 



where H p i VOt is the Hubble scale corresponding to fc p i V ot, 
which leaves the horizon 7V p i vot e-folds before the end of 



5 In practice, w is not exactly 1/3, if the number of degrees of 
freedom in the thermal bath is itself a function of temperature. 



6 The code also includes the capability to define the end of inflation 
as corresponding to a particular ifrend > which would be useful for 
implementing multi-field models. However, with the exception 
of the test in Sec. Ill El we do not make use of this feature in the 
present work. 



inflation. In what follows, the pivot scale is set to k p i vot = 
0.05 Mpc -1 . For a specific inflationary potential we can 
compute iVp ivot from the usual matching equation; the 
post-inflationary expansion is known by hypothesis once 
we assume instantaneous reheating and 



N$*> = 55.75 - log 



I0 16 GeV 

V 1/4 t 

pivot 



log 



v l/ \ 

pivot 



V 



1/4 



, (17) 



where the above expression is drawn from Eq. (20) of 
Ref. [6J], with appropriate substitutions. 7 In the IRH 

case we compute simply N^ ivot from this expression and 
the iVpivot is not an independent variable in the chains. 



In the GRH case we set the prior 20 < N pivot < N 



(IRH) 
Y pivot 



The lower limit comes from requiring that cosmologically 
relevant wavenumbers are far outside the horizon when 
inflation ends. This is tacitly assumed by our compu- 
tation of the spectrum from the asymptotic mode am- 
plitude in any case, and values of N near this limit are 
excluded by the data for all the models we consider here. 
The upper limit is enforced by rejecting any step to a 
model for which iV pivot > NJ^.' . 



E. Accuracy and timing 

To test the accuracy of ModeCode, we compare its 
output with the analytic solution for primordial pertur- 
bations in the "power law inflation" model [82| . where 
the scale factor evolves as a(t) oc t p during inflation and 
the potential has an exponential form 



V{4>) = V e*v 



pMp\ 



(18) 



Power law inflation is one of the very few known models 
for which the spectrum of primordial perturbations can 
be computed exactly, making it an ideal test case for the 
numerical solution of ModeCode. 

Since the tensor-to-scalar ratio for this model is r = 
16/p, we must take a large value of p to avoid violating 
present upper limits on r; here we choose p = 60. For 
p > 1, power law inflation does not end via slow roll 
violation [Eq. (|15|) ]. so we impose an end to the inflation- 
ary expansion at </> C nd = Mp\- We additionally assume 
log(Vb/Mpj) = —8.8 and -/V p i vot = 50, yielding power 
spectra that are reasonably consistent with observations. 

In Fig. []] we compare the ModeCode solution for 
scalar modes with the exact solution from Ref. [LI- At 
all scales, the numerical evolution of scalar perturbations 



7 This expression makes no allowance for the changing number of 
rclativistic degrees of freedom as the Universe cools or for dark 
energy, and it assumes a sharp transition between radiation and 
matter-dominated expansion. The resulting approximation will 
not significantly bias our results using present data, but may 
need to be addressed in the future. 
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FIG. 1: Test of the accuracy of ModeCode for power law 
inflation. Upper panel: comparison of the exact analytic so- 
lution (solid red curves) and the ModeCode solution (black 
points) for \uk\ as a function of conformal time r (following 
Ref. [13 |. r is negative during inflation). For comparison, the 
evolution of \z\ is plotted as a dashed green line. Lower panel: 
percent error in the ModeCode solution for each of the points 
plotted in the upper panel. 



matches the exact solution with an accuracy of about 
0.01% or better. This test indicates that ModeCode 
does not introduce significant error in the computation 
of CMB angular power spectra by CAMB, which has a 
root mean square accuracy of ~ 0.3% in the configura- 
tion to be used for the Planck analysis (A. Lewis, private 
communication) . 

In the course of previous work using ModeCode to 
compute numeric al p ower spectra for a potential with a 
step-like feature |59j, it was extensively compared with 
the independent code of Ref. [5(| , yielding agreement to 
numerical precision. It also agrees to similar precision 
with the mode evolution code used for the same potential 
in Rcfs. [53, [83[, once the nonstandard choice of initial 
conditions in the latter work is accounted for. 8 

To compute the primordial power spectra at arbitrary 
values of fc in CAMB, ModeCode uses cubic spline in- 
terpolation on a grid of k values spaced evenly in lnfc. 
The extra time required to run ModeCode with CAMB 
depends primarily on the number of k values in this 
grid. For the default setting of 500 grid points over 
10~ 5 < k/Mpc~ < 5, which provides more than suffi- 
cient accuracy for smooth primordial power spectra, us- 



See Appendix A in Ref. 15911 for details. 



ing ModeCode with CAMB typically requires < 15% 
more time per evaluation than the default version of 
CAMB. The number of grid points can be easily adjusted 
in the code to accurately deal with more complicated po- 
tentials for which finer sampling in k is required. 



F. Data 

CMB Data: We use the v4 version of the 7- yea r 
WMAP likelihood function with standard options 84|, 
the ACBAR bandpowers from Ref. (85[ between 550 < 
I < 1950, and the Pipeline 1 QUaD bandpowers be- 
tween 569 < i < 2026 from Ref. J86J]. For IRH mod- 
els, we consider constraints from the 7-year WMAP data 
("WMAP7") only. In the GRH case, we present both 
WMAP7 results and constraints that additionally include 
data from QUaD and ACBAR ("WMAP7+CMB"). 

Planck Simulation: For selected inflation models in the 
GRH case, we use an unpublished simulation kindly pro- 
vided by George Efstathiou and Steven Gratton. In this 
simulation, the model for the "observed" power spec- 
tra has four components: the primordial CMB power 
spectra, unresolved point sources, unresolved Sunyaev- 
Zel'dovich (SZ) clusters, and instrumental noise. The 
input CMB power spectra are computed from a random 
realization centered on the "best fit" WMAP 5-year cos- 
mological parameters, including n s = 0.963. In addi- 
tion the simulation includes a tensor component with 
r = 0.1, close to the margin of detcctability by Planck 
[871 l88j . The likelihood function is described by an exact 
Wishart distribution, marginalizing over the SZ model 
and the point source model as nuisance parameters in the 
MCMC. Note that while the fiducial model has r = 0.1, 
the particular realization used in the simulation is con- 
sistent with a somewhat larger tensor amplitude that is 
closer to r ~ 0.14. The best fit values for the models we 
consider reflect this larger tensor-to-scalar ratio. 



G. Models and priors 

In order to illustrate our methods, we derive con- 
straints on a variety of models. First we consider a se- 
quence of "single term" potentials, 



V = A — 
n 



(19) 



with n = 2/3, 1, 2, and 4. The last two cases correspond 
to the canonical quadratic and quartic chaotic inflation 
models [831, an d for n = 2 we replace A with m 2 in our 
discussion. String motivated scenarios [90, |91( can yield 
potentials with the form Eq. (|19[) and non-integer values 
of n at large </>. We assume that these potentials are mod- 
ified for 4> < to ensure that V(<fi) > 0, if n is not an even 
integer. For convenience, we will work with the simple 
monomial term; to explicitly constrain the corresponding 
stringy scenarios we would instead have to work with the 



TABLE I: Priors on model parameters and maximum likeli- 
hood (ML) values for WMAP7 GRH constraints. All GRH 
models include a uniform prior of 20 < N p i vot < N^ ivot . 
Dimensionful quantities are expressed in units where the re- 
duced Planck mass Mpi is set to unity. Values of n refer to 
specific cases of Eq. (119[) . 



Model 



Priors 



1s,ML fML 



-21n£ M L 



n = 2/3 
n= 1 
n = 2 
n — 4 
Natural 

Hilltop 



-11 < logA < -7.5 
-11 < log A < -7.5 
-12 <logm 2 < -8 



0.965 0.07 
0.969 0.08 
0.964 0.14 



-13.4 < log A < -10.4 0.949 0.27 
-5<logA<0 0.962 0.08 

0.5 <log/ <2.5 

-8<logA<-2.8 0.944 4x10" 
-13.3 < log A < -12 



7475.2 
7475.4 
7477.3 
7488.7 
7475.8 

7476.2 



full form of the potential, since the modification near the 
origin changes </> en d and the matching between k and <f>. 
For the single term potentials, the free parameter fixes 
the height of the potential and the amplitude of the per- 
turbations (i.e. A s ). However, for these models the other 
spectral parameters (n s , r, etc.) are well approximated 
by combinations of slow roll parameters, which do not de- 
pend directly on A. When the potential has two or more 
free parameters, the mapping between the explicit form 
of the potential and the cosmological observables grows 
more complicated, since these parameters affect both the 
height and shape of the potential. We give constraints 
on axion-motivated "natural inflation" K)2f with 



V(</>) = A^ 



+ cos [ - 



and "hilltop inflation" 93149a with 



A 



V^^A 4 --^ 



(20) 



(21) 



for which r takes almost arbitrary values while n s remains 
close to unity. 

We show the model priors used in our MCMC analysis 
in Table [H The model parameters correspond to un- 
known scales in high energy particle physics, so it is nat- 
ural to sample them logarithmically. For the single pa- 
rameter models, we will see that the data constrain these 
scenarios more strongly than the priors on the height 
of the potential. On the other hand, the priors chosen 
for the natural and hilltop inflation models are more re- 
strictive. Both natural and hilltop inflation have limits 
in which they are essentially identical to a <p n model, 
and one of their two free V(<j)) parameters is irrelevant 
[6J, uM- For natural inflation, this limit is / — s- oo, such 
that 
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FIG. 2: Constraints on the quadratic potential (n = 2). 
Left: constraints on logm 2 and N p i vot for WMAP7 (large 
gray contours) showing 68% CL (light shading) and 95% CL 
(dark shading) contours. Green dashed curves show the con- 
tours for WMAP7+CMB data, and the small blue contours 
show simulated Planck constraints. Right: marginalized ID 
n s and r distributions for WMAP7 (thick solid curve, black), 
WMAP7+CMB (dashed curve, green), and the Planck simu- 
lation (thin solid curve, blue). In both panels, the sharp right 
edges of the distributions correspond to the IRH constraints, 
for which the values of iVpi vo t (left), n a , and r (right) are 
nearly fixed. Here and in Figs. [5] and [6] the input parameters 
for the Planck simulation are n s — 0.963 and r — 0.1 (see 
Sec.HTFl). 



FIG. 3: Same as Fig. [2]for the quartic potential (n = 4). 



with 0o = n fy which is a quadratic potential after a 
field redefinition. In the hilltop case, if A is very large, 
the astrophysically relevant portion of the potential is far 
from the origin. In this limit, 



V((t>) « 4A 4 1 - f- 



(23) 



where O = v^A/A 1 / 4 is the field value at which the po- 
tential crosses zero. A field redefinition yields a purely 
linear potential with a single free parameter. Our pri- 
ors are chosen to avoid the regions of parameter space 
where the two parameters of the natural or hilltop infla- 
tion model are degenerate. 

We stress that the goal of the current paper is to in- 
troduce ModeCode, explore the ability of both current 
and anticipated datasets to distinguish between differ- 
ent inflationary models, and constrain the number of e- 
folds of inflation required by specific inflationary scenar- 
ios. Broader issues surrounding the choice of priors and 
model selection will be addressed in a forthcoming paper. 



III. RESULTS 

A. Single term potentials 

Figures EHU show the constraints from WMAP7 and 
WMAP7+CMB on log A (or logm 2 for n = 2) and N pivot 
for V(4>) = \<f> n /n with n = 2/3, 1, 2 and 4, along with 
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FIG. 4: Constraints on single term potentials with n = 1 
(top) and n — 2/3 (bottom). Conventions match Fig. [21 



the inferred constraints on n s and r. The limits on log A 
from current data are much stronger than the priors listed 
in Table |U Stipulating instant reheating (IRH) ensures 
that A^pivot is nearly independent of the other model pa- 
rameters, and only weakly dependent on the value of the 
exponent n: N v - lvot ' ~ 57-59. In the general reheating 
(GRH) case, slow roll calculations lead us to expect that 
A^pivot, n 3 , and r are all effectively functions of a single 
free parameter, and are thus strongly correlated; these 
correlations arc reflected in the more accurate Mode- 
Code constraints. The IRH upper limit from the prior 



on K 



pivot 



creates a sharp cutoff in the distributions for 



n s and r. Given the prior on 7Vpi vot , each of these mod- 
els has a red tilt (n s < 1) and a "large" tensor-to-scalar 
ratio, r ~ 0.1. 

If ATpivot is lower than its IRH value, or w < 1/3, A (or 
m 2 ) is larger than its IRH value: the overall range in this 
parameter is typically a factor of ~ 2-3. Lower values of 
-^Vpivot correspond to smaller n s and larger r, as does 
increasing the exponent n for fixed iV p i vot . Specifically, 
the slow roll approximation gives 



N, 



pivot 



C 



— n — 

2 n-1 



(1-n,)- 1 ^— , (24) 
r 



where £ = for n = 1 and ( = 1 otherwise. Thus upper 
bounds on 1 — n s and r from data can set lower bounds 
on A^pjvot. Except for the n = 4 case, the predicted value 
of r is less than the current WMAP upper bound for 
all values of iVpi vot allowed by the prior. For n = 2/3 
and n = 1 the constraint on N p i vo t appears to be driven 
largely by the correlation between this parameter and 
n s — in all cases n s < 0.94 is strongly excluded. For 
n = 2, the larger value of r found with smaller Af p i vo t 
provides some additional constraining power. Including 
CMB data on smaller angular scales from QUaD and 
ACBAR slightly strengthens these limits, relative to the 
WMAP 7 constraints, but does not significantly alter our 
conclusions. 

Of the single term potentials we consider, the quartic 
\<t> A /A potential has the largest tensor amplitude and the 
greatest deviation from scale invariance. In agreement 
with previous analyses of CMB data (e.g. Ref. |5l|). 
this model is excluded by WMAP7 data. Specifically, 
— 21ii£ml exceeds the values found for all of the other 
models considered here by > 12 (see Table [XJ> . Superfi- 
cially, this does not appear to be significantly stronger 
than the result obtained with a single year of WMAP 
data [51J. However, the analysis of Ref. [51j was carried 
out at fixed Af p i vo t = 50, whereas Af p i VO t is a free pa- 
rameter in our chains. If we impose the additional prior 
A^pivot < 50, the maximum likelihood of the quadratic 
potential worsens to — 21ii£ml = 7499.4. This is ~ 24 
larger than the overall best fit and excludes such mod- 
els with much greater confidence than the first year of 
WMAP data alone. 

In recent work, Martin and Ringeval [631 ] quote con- 
straints on the reheating temperature following inflation 
driven by a V((j>) ~ <p n potential. The tightest con- 
straints they present imply that the reheat temperature 
is above the TeV scale. However, this specific constraint 
is obtained for a prior that renders the post-inflationary 
expansion rate a function of n. Since scenarios for which 
V(4>) ~ </>" at large field values can have very different 
shapes near the origin, the prior could only be realized 
by a carefully tuned potential, as V(<j)) would need to 
be well approximated by <p n at energies far below the in- 
flationary scale. Moreover, this form of V(4>) must be 
modified near the origin if V(4>) > for all <fi, and the 
potential does not have a discontinuous first derivative 
at <f> — for n ^ 2 or 4. 
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FIG. 5: Top left: WMAP7 68% and 95% CL constraints 
on natural inflation parameters log A and log / for the GRH 
(large contours, gray shading) and IRH (small contours, red 
shading) scenarios. Top right: constraints on n s and r; points 
show random samples of models from the MCMC analysis 
for which -/V p i vo t is within 0.25 of the values indicated in 
the plot. Bottom: Predicted natural inflation GRH con- 
straints from Planck (small contours, blue shading) compared 
with current constraints from WMAP7 (gray shading) and 
WMAP7+CMB (dark green, dashed contours). 



Constraints using the simulated Planck likelihood show 
that the uncertainties in parameters of these single term 
potentials will be greatly reduced by the next generation 
of cosmological datasets. For the quadratic potential, 
we see from Fig. [2] that the simulated Planck constraints 
strongly disfavor (at > 95% CL) models with n a < 0.96, 
corresponding to N p i vot < 50 and logm 2 > —10.3. Re- 
call that for quadratic inflation followed by a matter- 
dominated phase and then thermal inflation, N p i vot is 

at least 10 less than the instant reheating value N^ iyot 
[761 ]. Consequently, we predict that Planck can differ- 
entiate between these two post-inflationary scenarios for 
quadratic inflation. 



B. Natural inflation 

Figure [5] shows our constraints on the natural inflation 
[Eq. (|2"0")) ] parameter space and the derived empirical pa- 
rameters n s and r from ModeCode. The relationship 
between the empirical parameters and the potential pa- 
rameters for natural inflation is discussed in detail in Ref. 
[961 ] . along with parameter constraints derived from the 



10 



3-year WMAP dataset. 

Unlike the single term potentials, current data permit a 
wide range of natural inflation parameters and, as noted 
in Sec. Ill Gl there is a degeneracy between / and A in 
the limit where these parameters are large. In this re- 
gion of parameter space, natural inflation overlaps with 
the quadratic model. Our priors are chosen to exclude 
most of this region; given the parametrization of the po- 
tential, if we allowed arbitrarily large values of / and A 
(and given that the quadratic potential is not currently 
excluded by data) almost all points drawn by the chains 
would be in this degenerate region. Our adopted priors 
on log / and log A still allow a region of nearly-degenerate 
models that contribute to the "ridge" seen in the right 
panels of Fig. [5] these models closely match the values of 
n s and r seen in the quadratic potential constraints. The 
marginalized constraints on n s and r depend strongly on 
the prior on log / due to the projection of a large number 
of degenerate models into this ridge. Thus, the apparent 
preference for this region of parameter space over models 
with lower values of r is largely due to this effect, and is 
not driven by the data. 

Instant reheating requires A pivot ' ~ 58, similar to the 
constraint for the single term potentials, but the addi- 
tional inflationary degree of freedom in the potential per- 
mits a larger range of n s and r. More generally, for fixed 
-^pivotj decreasing A and / reduces both n s and r. This 
is shown by the MCMC samples with fixed A p i VO t plotted 
in the upper right panel of Fig. [5] Thus natural inflation 
models can have lower values of r than the quadratic 
potential in Fig. [21 without increasing n s and N p i vot . 

The lower panels of Fig. [5] show how the natural in- 
flation constraints improve with additional CMB data. 
As for the single term potentials, the difference between 
WMAP 7 and WMAP 7+ CMB constraints is small, but 
Planck is expected to yield a dramatic improvement. In 
particular, the uncertainty in A p ; vo t — which is visible in 
the width of the log A — log / contours — is substantially 
reduced in the Planck forecast. Since the Planck simula- 
tion we use is not inconsistent with m 2 <f> 2 inflation, the 
limit in which natural inflation becomes indistinguishable 
from quadratic inflation would not be excluded in this 
particular forecast. Conversely, if quadratic inflation is 
disfavored by future data, we will be able to put data- 
driven upper bounds on / and A in addition to tightening 
the existing lower bound. 



C. Hilltop inflation 

Figure [6] shows the constraints on log A, log A, n s , and 
log r for hilltop inflation. We impose an upper limit 
of A < 0.0015 to remove models where the field starts 
far from = and near the V — crossing point, 
as described in Sec. Ill CI The remaining models have 
a small tensor amplitude and relatively large deviations 
from scale invariance. The MCMC samples plotted at 
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FIG. 6: Same as Fig. [5] for the hilltop inflation model with 
parameters log A and log A. Note the logarithmic scale for r 
in the right panels; the prior here permits very small values 
of r. 



correlated with n s , but not r, over most of the allowed 
region of parameter space. 

Constraints on the hilltop inflation model from 
WMAP7+CMB data and from the Planck simulation are 
compared with the (GRH) WMAP7 constraints in the 
lower panels of Fig. [6j The values of n s and r allowed 
by the hilltop model within our chosen priors are smaller 
than those assumed in the Planck simulation, so the con- 
tours for Planck are concentrated at the largest allowed 
values of both parameters. In fact, for this particular 
forecast the entire region of the hilltop inflation param- 
eter space within our priors would be strongly excluded 
by Planck, with — 2\tiCml ~ 75 larger than its value for 
the quadratic and natural inflation models. 

Unlike the other models considered here, current data 
allow jVp ivot ^ to cover a substantial range (roughly 10 
e-folds), as illustrated in Fig. [7] The weak constraint 
on A pivot / is due to a special cancellation in the slow 
roll expression for the scalar spectral amplitude, which 
leaves A s independent of the overall height of the po- 
tential. 9 Meanwhile, the departure from scale invariance 
is 1 — n s ss 3/Ap ivot in the limit of small A. Other 
two-parameter models can formally support inflation at 
low energy scales, but the spectra of such models are 



fixed Apivot in Fig. [5] show that the number of e-folds is not generic. 



9 In the more general class of hilltop potentials V{4>) = A 4 — \<f> n /n, 
this situation only occurs for n = 4 [93H95I ]. so this property is 
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FIG. 8: The parameter space in n 3 and r corresponding to 
uniform sampling within the priors on potential parameters 
listed in Table [T] (unshaded curves, 68% and 95% CL regions) 
for natural inflation (left) and hilltop inflation (right). The 
WMAP7 GRH constraints from Figs.[5]and[6]are shown again 
for comparison as shaded gray 68% and 95% CL regions. Or- 
ange shading shows the 95% CL region for empirical WMAP7 
constraints with flat priors on n s and r. 



typically far from scale-invariant and thus disfavored by 
the data. On the other hand, hilltop inflation can have 
A ^C 10 16 GeV, which ensures that the tensor amplitude 
is very low, without driving the spectral index to an ob- 
servationally excluded value. Consequently, both r and 



N. 



(IRH) 



pivot 



can vary greatly, as we see in Fig. [6] Models with 
very low values of A do have a lower likelihood, as the 
stronger breaking of scale invariance in these models is at 
odds with the spectral tilt allowed by the data. For IRH 
models there is effectively only one free parameter in the 
potential after fixing A s , leading to a strong correlation 
between n, and r which is absent in the GRH case. 



D. Impact of model priors 

For the single term potentials, the priors on the ampli- 
tude of the potential A (or m 2 ) listed in Table U are suffi- 
ciently weak compared to the constraints from the data 



that they have no significant impact on the estimated pa- 
rameter values and confidence regions. The upper limit 
on iVpi VO t corresponding to instant reheating, however, 
does significantly reduce the allowed region in parame- 
ter space for each of these models, thus restricting the 
possible values of n s and r. 

The resulting upper limit on n s and lower limit on 
r can lead to tension with observations. For example, 
the WMAP7 constraints are n s = 0.982t° 2 1 9 (68% CL) 
and r < 0.36 (95% CL) [84||, treating these as empiri- 
cal parameters without specifying a particular potential. 
For the quartic A0 4 /4 potential, the instant reheating 
limit on iV p i vo t and the nearly perfect correlation be- 
tween A^pivot, n s , and r result in the limits n s < 0.95 
and r > 0.27; since both of these are in tension with the 
measured values, the quartic potential provides a poor 
fit to the data relative to the other single term potentials 
which can achieve larger n s and smaller r. 

On the other hand, as noted earlier the natural infla- 
tion and hilltop inflation priors have been chosen specif- 
ically to limit the extent of parameter degeneracies that 
would otherwise be allowed by current data. In the lim- 
iting regions that are truncated by the priors, these two 
models are degenerate with specific monomial potentials, 
as discussed previously. Furthermore, the mapping be- 
tween 8y and 6 cmp for these models is significantly more 
complicated than for the single term potentials, so the 
uniform top-hat priors on 0y can correspond to highly 
non-uniform prior distributions for 6 cnlp . This mapping 
of the priors is illustrated in Fig. which shows the re- 
gions of the (n s , r) plane obtained by uniform sampling 
within the priors on 6y specified in Table[T]for the natural 
inflation and hilltop inflation models. 

Due to the natural inflation degeneracy between log A 
and log/, even in the absence of any data the priors 
clearly favor models along the line r w 4(1 — n s ) cor- 
responding to the approximately quadratic regime near 
the minimum of the natural inflation potential. The 
WMAP7 GRH constraints in Fig. [5] are qualitatively de- 
scribed by the intersection of the priors with the empiri- 
cal constraints on n s and r from WMAP7. Note that the 
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upper limit on n s is typically set by the -/V p i vot < N pivot ' 
prior (the lower right edge of the 95% CL region in the 
left panel of Fig. [5]) , while the lower limit on n s — and 
thus iVpi vot — follows from the data. The 95% lower limit 
on r does not quite reach down to the limit allowed by 
the prior because of marginalization over the strong pro- 
jection effect described in Sec. lIIIBl which favors models 
along the upper diagonal ridge of the prior. With present 
data, the region in the (n s ,r) plane allowed by the pri- 
ors is fully consistent with the measured values of these 
parameters (and therefore the best fit n s and r values in 
Table Q] arc consistent with the WMAP7 empirical con- 
straints), but upcoming measurements may yet rule out 
the natural inflation model. 

For the hilltop inflation model, the constraints on n s 
and r are partially influenced by the priors on log A and 
log A from Table HJ in particular, the range of r allowed 
is limited by the prior on log A (see Fig. [5J. However, 
the distribution of models allowed by the priors in the 
(n s , log r) plane is much more uniform for hilltop inflation 
than it is for natural inflation. As for natural inflation, 
the upper limit on n a for hilltop inflation models is set 
by the A p i VOt < N pivot ' prior, and the lower limit on n s 
is determined by the data. Like the quartic potential, 
the upper limit on n s corresponding to instant reheating 
is low compared with the preferred value from WMAP7; 
however, in the case of hilltop inflation, this is coupled 
with a small value of r, which enables the hilltop inflation 
model to fit the WMAP7 data reasonably well since the 
empirical constraints on n s and r are correlated. 
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FIG. 9: Approximate constraints on the quadratic potential 
obtained by selecting MCMC samples from standard WMAP7 
ACDM+tensor chains and mapping A s , n s , and r to login 2 
and JVpivot using the slow roll approximation. The shaded 
gray contours are the WMAP7 GRH constraints from Fig. [2] 
obtained by a direct MCMC estimate of m 2 and A p i vo t for 
quadratic inflation. The unshaded orange contours are de- 
rived by sampling chains for A s , n 3 , and r, from which we 
selected only those points with \1 — n s — r/4| < 0.005 to ap- 
proximately match the slow roll relation between n a and r for 
the quadratic potential. 



E. Slow roll mapping 

Many previous constraints on inflationary models have 
been obtained by taking the empirical parameters cmp 
(e.g. A s , n s , and r) predicted by a given model and com- 
paring these with constraints on empirical parameters 
(see e.g. [ll,[52|)- Given current data, limits on A s , n s and 
r typically dominate the constraints on simple inflation- 
ary models; starting from constraints on these empirical 
parameters allows constraints on a number of inflation- 
ary models to be inferred from a single set of empirical 
Markov chains. 

For the simplest inflationary models, this approach al- 
lows one to quickly constrain many potentials at once. In 
general, however, the mapping B cmp — > By is not one-to- 
one; thus this method requires sampling within the space 
of potential parameters and is no more efficient than di- 
rectly constraining these parameters from the data using 
ModeCode. For any potential parametrized by By, one 
can compute the primordial scalar and tensor power spec- 
tra (as described in Sec. Ill Al or otherwise) and find the 
values of 8 omp by fitting to these power spectra. However, 
for any given inflationary potential many values of 6 cmp 
cannot be obtained for any combination of By . Likewise, 
as we saw for the natural and hilltop scenarios, some 
models have degenerate combinations of By which cor- 



respond to the same primordial power spectra, so the 
mapping By — > 6 cmp is not always invertible. Further, 
models with sharp features or other complexities do not 
produce a power spectrum that is easily described by the 
usual empirical paramctrization. 

Despite these problems, it is instructive to attempt to 
derive constraints on By from constraints on B euip for 
the purposes of comparison with our main results in the 
previous sections. Here we perform these tests using the 
slow roll approximation, although one could implement 
the mapping using more precise methods. In Appendix |A1 
we describe the slow roll mapping for the quadratic po- 
tential. Figurc[9]shows that the constraints on log m 2 and 
A^pivot obtained by mapping from empirical parameters 
agree reasonably well with the direct MCMC constraints 
from ModeCode. However, they are visibly noisy due 
to the smaller number of MCMC samples, and are sys- 
tematically shifted toward lower values of log m 2 , due to 
the mapping being done at lowest order in slow roll. 

While this approach works well for the quadratic po- 
tential (and can be expected to produce similar results 
for the n = 1 and n = 2/3 cases) it is far less efficient 
for the quartic (n = 4) potential. Recall that this sce- 
nario is disfavored by the data. Consequently, using the 
Metropolis-Hastings algorithm to draw samples from the 
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(n s ,r) parameter space results in very few (if any) ac- 
cepted points in the relevant region for the quartic po- 
tential on the (n Sl r) plane, and thus the contours of A 
and ./Vpivot are extremely noisy. Constraining the quartic 
potential parameters using this method requires MCMC 
sampling that is specifically designed to acquire samples 
in the (n s ,r) region spanned by this model. 

Models with multiple parameters like natural inflation 
and hilltop inflation face an additional problem in that 
the mapping from By to O mp is not invertible in regions 
where there is a parameter degeneracy. For example, at 
fixed values of A 4 / -2 , all natural inflation models with 
large / have nearly identical power spectra (see Sec. lIIGI) . 
Therefore, a single point in 9 cmp space cannot be simply 
mapped to the corresponding values of By in this degen- 
erate region. Such a mapping would require an additional 
MCMC run or some similar method of sampling in the 
By parameter space. 



IV. DISCUSSION 

In this paper we introduce ModeCode, a new, pub- 
licly available numerical solver for the inflationary per- 
turbation equations. ModeCode is extendable, numer- 
ically efficient, and integrated with CAMB and Cos- 
moMC. We demonstrate the use of ModeCode by con- 
straining several single-field inflationary models using 
current CMB data, confirming that present-day data put 
useful constraints on the parameters of a number of in- 
teresting and well-motivated inflationary models. 

Using a simulated likelihood, we present forecasts for 
the quality of the constraints that can be expected from 
Planck, showing that these will greatly strengthen limits 
on the inflationary parameter space and possibly exclude 
some simple inflationary models. In the Planck forecast 
analysis, our aim was not to provide "Fisher" -style fore- 
casts for each model in turn, but to take the realiza- 
tion provided by the simulation and analyze it with our 
pipeline as we would the real sky. Given the level of ten- 
sors present in this simulation, which is at the margin of 
detectability with Planck, we find that small-field mod- 
els would be excluded with high significance as expected. 
Conversely, in a scenario where Planck polarization mea- 
surements did not find evidence of tensor fluctuations, 
we predict that many large-field models would be cither 
excluded or limited to a small region of parameter space. 

In our analysis we pay close attention to the inter- 
play between the post-inflationary expansion history and 
the inflationary observablcs. These arc connected by the 
matching criterion, which determines the moment dur- 
ing the inflationary epoch at which a given comoving 
scale leaves the horizon. Inflation is often assumed to be 
a GUT scale phenomenon, but the expansion rate and 
thermal state of the post-inflationary universe is not di- 
rectly constrained until McV scales, at which point the 
success of Big Bang nucleosynthesis and evidence for a 
cosmological neutrino background strongly suggest that 



the universe was thcrmalizcd. Consequently, there is a 
huge range of energies over which the composition and 
expansion rate of the universe are effectively undeter- 
mined. 

For most inflationary models, the constraints on the 
total number of e-folds since the pivot scale left the hori- 
zon, A^pi vo t, are noticeably tighter than the hard lower 
bound (A/pivot > 20) assumed in our analysis. Thus we 
can be confident that our constraints on iVpivot ar c driven 
by the data, although these currently eliminate only rel- 
atively extreme post-inflationary scenarios. By contrast, 
the Planck forecast suggests that the next generation of 
CMB data will put much tighter constraints on the re- 
heating history, given a specific inflationary model. For 
instance, for m 2 4> 2 inflation, Planck should discriminate 
between a long matter-dominated phase that extends to 
the TeV-scale and instant reheating at perhaps the 2 a 
level. 

This correlation between the post-inflationary dynam- 
ics and the inflationary epoch has significant conse- 
quences for particle physics. For instance, many super- 
symmetric scenarios predict that the primordial universe 
undergoes a period of matter domination driven by heavy 
moduli (e.g. [78|, |97|). for which A p i vot differs substan- 
tially from the instant reheating value, and post-Planck 
cosmology will thus be increasingly concerned with the 
full evolutionary history of the universe. Moreover, infla- 
tionary model builders will be able to definitively test sce- 
narios which include predictions for the post-inflationary 
expansion rate. 
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Appendix A: Slow roll mapping — quadratic potential 

Here we describe the procedure that produced the con- 
straints on the quadratic potential shown in Fig. |H] by 
using slow roll relations to map constraints on the em- 
pirical parameters 9 cmp = {h\A s , n s , r} to constraints on 
the potential parameters By = {A p i vot ,logm 2 }. For the 
potential V = m 2 <p 2 /2, the slow roll approximation gives 



m 

-^Vpivot 



24n 2 A s e 2 v , 
-A^-n s ) 



16' 



(Al) 
(A2) 

(A3) 



where e v = {M^/2)(V^/V) 2 . 

We begin with constraints on 6 cmp from a standard 
MCMC analysis of WMAP7 data, which treats A s , n s , 
and r as free, independent parameters. We then post- 
process the empirical parameter chains to express the 
constraints in terms of 8y by first imposing new priors 
on 8 cmp , and then using Eqs. (|A1|) — (|A3|) to compute the 
values of 9y for each MCMC sample. 

The first prior applied to C mp enforces the slow roll 
relation between n s and r for the quadratic potential, 
r m 4(1 — n s ), by selecting only those points from the 
chains that satisfy |1 — n s — r/4| < 6. The parameter 
(5 can be adjusted either to match the slow roll relation 
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more closely in the limit S — > 0, or to retain more MCMC 
points for better sampling of the likelihood function (and 
thus smoother contours) by choosing a larger value of 
5; for Fig. [5] we have chosen S = 0.005. Note that the 
allowed values of A s are not restricted by specializing to 
the case of the quadratic potential. 

We additionally exclude points for which iVpi vot is large 
enough that prh is greater than (3/2)Vend, thus violating 
energy conservation. (The factor of 3/2 arises because 
the energy density at the end of inflation contains a sig- 
nificant contribution from the kinetic energy.) That is, 



we require iVj 



pivot 



<JV. 



(IRH) 



pivot 



where 



A, 



(IRH) 
pivot 



ln(H P ivot/fcpivot) - 71.1 - M^/MpO, 



(A4) 



with Vend ~ Jn 2 Mpj for the quadratic potential. 

For the MCMC samples that remain after applying 
the cuts on n s , r, and A p i vot , we multiply the likelihood 
of each sample by r~ 2 to approximate flat priors on 8y 
instead of the original flat priors on cmp . 

After including each of these priors and mapping the 
MCMC parameters from 8 cmp to 9y using Eqs. (|A1I) 
(|A3[) , the resulting constraints shown in Fig.[5]match the 
more accurate constraints from ModeCode (Sec. IIII A[) 
reasonably well, albeit with noisy contours due to poor 
sampling along the line r ~ 4(1 — n s ). The remaining 
systematic offset between the two sets of constraints is 
consistent with the effects of second order corrections to 
the slow roll approximation. 
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